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ABSTRACT 

The physical mechanisms and timcscales that determine the morphological signatures and the 
quenching of star formation of typical (~ L*) elliptical galaxies are not well understood. To ad- 
dress this issue, we have simulated the formation of a group of galaxies with sufficient resolution to 
track the evolution of gas and stars inside about a dozen galaxy group members over cosmic history. 
Galaxy groups, which harbor many elliptical galaxies in the universe, are a particularly promising 
environment to investigate morphological transformation and star formation quenching, due to their 
high galaxy density, their relatively low velocity dispersion, and the presence of a hot intragroup 
medium. Our simulation reproduces galaxies with different Hubble morphologies and, consequently, 
enables us to study when and where the morphological transformation of galaxies takes place. The 
simulation does not include feedback from active galactic nuclei showing that it is not an essential 
ingredient for producing quiescent, red elliptical galaxies in galaxy groups. Ellipticals form, as sus- 
pected, through galaxy mergers. In contrast with what has often been speculated, however, these 
mergers occur at z > 1, before the merging progenitors enter the virial radius of the group and before 
the group is fully assembled. The simulation also shows that quenching of star formation in the still 
star- forming elliptical galaxies lags behind their morphological transformation, but, once started, is 
taking less than a billion years to complete. As long envisaged the star formation quenching happens 
as the galaxies approach and enter the finally assembled group, due to quenching of gas accretion and 
(to a lesser degree) stripping. A similar sort is followed by unmerged, disk galaxies, which, as they 
join the group, are turned into the red-and-dead disks that abound in these environments. 
Subject headings: galaxies: elliptical and lenticular, cD - galaxies: evolution - galaxies: groups: general 
- galaxies: interactions 



1. INTRODUCTION 

Galaxy groups in the local universe harbor a large 
population of massive galaxies with elliptical morpho l- 
ogy and quenched star formation (jKovac et al.1 l2010h . 
Observations indicate that, to a large extent, this re- 
sults from environmental forcing o ver the past eigh t 
billion years since redshift z ~ 1 (|Peng et all 12010( 1 . 
Many physical processes may be responsible for mak- 
ing elliptical morphologies and quenching star forma- 
tion; major and minor merging, ram pressure strip- 
ping, accretion shocks, and removal of hot gas are 
among those that have b e en suggested as possible 
mechanisms (iGunn fe Gottl Il972t iLarson et all 119801: 
QuilisetalJ 120001: iBalogh fe Morris! 120001: lCox_et_al 
20061: iBirnboim fc Dekell |2003t Kawata fe Mulchaev 



200a iMcCarthv et all l200alHopkins et al.ll2009f l. The 
challenge is to connect the large-scale, cosmological for- 
mation of the group, which affects the global evolution of 
group members, to the galactic sub-kpc scale, where star 
formation and stellar feedback do minate the physics of 
the gas and the galaxy morphology ((Springel et alll2"o"05t 
IGovernato et alj|2007l) . 

Both numerical simulation s of the galaxy population in 
a cosmological context(e.g., iKawata fc Mulchaev! 120081 : 
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iSimha et al.ll2009t ICroft et al.ll2009l ) and semi-analytical 
models of galaxy formatio n (e.g., iBenson fc Devereuxl 
l2010t iHopkins et al.l l2010al rbT) have, however, lacked so 
far the ability to describe at once the various environ- 
mental processes that can affect the baryonic compo- 
nent of galaxies, due to lack of resolution or the incom- 
plete knowledge of the interplay between such processes. 
Therefore, despite the observational evidence that galaxy 
evolution is affected by the galactic environment, how 
this happens within the cold dark matter scenario for 
the formation of cosmic structure is highly debated to- 
day. This notwithstanding the results of detailed, small 
scale numerical models, ab stracted from the large-scale 
cosm ological context (e.g., lOuilis et al.ll2000t ICox et all 
2006), have repeatedly shown how mergers, tidal inter- 
actions, ram pressure stripping of the cold and hot gas 
phase as well as gas accretion from the surrounding inter- 
galactic medium can change the morphology of individ- 
ual galaxies, and alter their star formation rates (SFRs) 
and colors. 

Recent progress in the modeling of the formation of 
galaxies in the cosmological context has finally shown 
that realistic galaxies can be formed in the cold dark mat- 
ter model once galaxies are resolved wit h ~ 10 5 baryonic 
particles and at a sub-kpc resolution (IGovernato et all 
l2TMlN^a^etld1l2W7l:lGuedes et alll201lH . but this has 
only been achieved so far for individual galaxies in low- 
density (field) environments. 

Therefore, key questions remain when, and how, the 
disks are transformed into spheroids, star formation is 
quenched, and stellar mass is accreted in the massive 
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galaxies that end up residing preferentially in the dense 
environments at z = 0. We address these key ques- 
tions with the help of a cosmological hydrodynamical 
simulation that follows the concurrent evolution of many 
galaxies before, during, and after they join a relatively 
dense, group-sized virialized region with virial radius 
-R200 = 396 kpc and virial mass M200 = 1.1 x 10 13 M 
at z ~ 0. 

The structure of this paper is as follows. We de- 
scribe our simulation and analysis methods in Section [5] 
and [3l respectively. In Section [4] we present and discuss 
our main results, including the properties of the group 
satellite galaxies (Section l4.ll and I4.2|) . and the origin 
of the morphological and photometric transformations 
(Section I4.3|) . We conclude in Section [5] 

2. SIMULATION DETAILS 

The simulation set-up and the employ ed zoom-in tech- 
nique is described in detail elsewhere ()Feldmann et al.l 
2010). In brief, a group halo is selected from a 
dark matter simulation of a 123 Mpc periodic box. 
Higher resolution initial conditions centere d on the group 
with WMAP-3 cosmological param eters (ISpergel et al.1 
|2007|) are generated using graf ic-2 (jBertschingerl l 200 1 ) . 
The simul ation is run with th e parallel TreeSPH code 
GASOLINE (jWadslev et al.ll2004D . The mass (force) reso- 
lution of the simulation is 6.4 x 10 6 M (0.5 kpc), 1.4xl0 6 
M Q (0.3 kpc), and 4 x 10 5 M Q (0.3 kpc) for dark matter, 
gas and star particles, respectively, which allows us to re- 
solve the half-mass radii of the massive group members 
(~ 1 kpc). We output 92 snapshots between z = 41.5 
and z = 0.1 equidistant in time, i.e., with At ss 130 Myr. 

Besides following gravitational, hydrodynamical, and 
optically thin radiative processes of matter in an expand- 
ing universe directly, the simulation includes modeling 
of star formation, feedback from type la and type II su- 
pernovae, mass loss by stellar winds, and metal enrich- 
ment. The simulation parameters were or iginally tuned 
to produce disk galaxies in g alactic halos (|Stinson et all 
l2006t iGovernato et al.ll2007t ). We use a star formation 
efficiency per free fall time of 0.05, and allow stars to 
form in a probabilistic fashion if the density exceeds 0.1 
cm~ 3 , the gas temperature is lower than 15,000 K, the 
gas is in an overdense (S > 55) region, and part of a 
convergent flow. A thermal feedback of 4 x 10 50 erg 
is injected into the gas per supernova. The cooling of 
the gas is turned off for the time corresponding to the 
end o f the snowplow phas e of a supernova type II blast 
wave ( Stinson et al.l [2006h . The cooling is not disabled 
for supernovae of type la. The simulation includes stel- 
lar mass loss due to stellar winds and a sp atially uniform 
UV background ()Haardt fc Madaul [1996V The simula- 
tion we present here c orresponds to the run G2-HR in 
iFeldmann et al.l (|2010D where it has been used to ana- 
lyze the properties of a massive, central group galaxy. 

As we show in Section 14.21 the same physical prescrip- 
tions produce, when used to simulate a relatively massive 
potential of a group of galaxies, z ~ galaxies with a 
broad range of luminosities, sizes and kinematics, which 
are roughly consistent with the properties of disk and 
elliptical galaxies in the local universe. They are also 
adequate to produce mas sive, early- type galaxi es at the 
centers of galaxy groups ()Feldmann et al.ll2010f) . 



3. HALO EXTRACTION AND THE MEASUREMENT OF 
THEIR PROPERTIES 

Halos and halo c enters are identified with the AMIGA 
HaloFinder (AHF; iGill et al.|[200l iKnollmann fc Knebd 
2009), that is capable of detecting sub-halos (our satel- 
lites) within halos (our group). Halos extend out to the 
virial radius -R200, which we define as the radius enclosing 
a mean matter density of 200/57 m times the mean density 
of the universe, or out to the tidal radius, whichever is 
smaller. The virial mass M200 is the mass within i?200- 
We then use the following technique to track satellites 
over time. First, we identify each halo h in each snapshot 
s below z = 4 that (1) has a bound mass of more than 
1.4 x 10 10 M Q , (2) contains more than 50000 particles, 
and (3) lies within 630 comoving kpc from the group cen- 
ter. Second, for each such determined halo h s we deter- 
mine its main progenitor h s -% in the previous snapshot 
and its main successor h s+ \ in the following snapshot and 
construct a complete timeline (ho, ■ ■ ■ , h s , . . . , h n ) by re- 
cursion. Let Ni, N2 and A common be the number of parti- 
cles in halo hi, the number of particles in halo /12 and the 
number of particles common to both halos, respectively. 
We then use the ratio 77(^1, /12) = N^ ommon / (NiN 2 ) 
to define main successors and progenitors. Specifically, 
we define the halo with the largest 7/ from snapshot 
s + 1 (any halo is in principle eligible, not only the 
ones selected as described above) as the main successor 
of halo h s and the halo with the largest 77 from snap- 
shot s — 1 as the main progenitor (but excluding the 
parent group halo). Since the group halo ho is a po- 
tential main successor this algorithm sometimes finds 
multiple timelines for a given halo h s . For instance, 
the following timelines of halo h s could be produced: 
(ho, ■ ■ ■ , h s , h s+ i, . . . , h n ), (ho, ...,h s , h s +i,hG, ■ • • , ho), 
and (ho, . . . ,h s , he ■ ■ ■ , ho). In this case, we only keep 
the timeline that traces the halo h s as long as possible 
as an entity separate from the group halo (which is the 
first timeline in the example). Following this procedure, 
we obtain ~ 20 satellites of which 13 have not merged 
with the central galaxy by z = 0.1. 

Unless noted otherwise SFRs are measured within a 
sphere of 20 physical kpc. The specific SFR is obtained 
by dividing the SFR by the stellar mass. Stellar masses 
are defined as the mass of the stellar component that is 
bound to the halo (using AHF) within a sphere of 20 
physical kpc. This radius encloses almost all (~ 98% on 
average, 91% for the largest galaxy at z = 0.1; ~ 100% at 
z > 1) of the stellar mass of non-central group members 
within their virial (or truncation) radius. 

Kinematic properties are determined from a slit that 
is 24 kpc long, 2 kpc wide and contains 24 bins along 
the two-dimensional major axis of an edge-on view of 
the galaxy. The rotation velocities and dispersions are 
determined from the moments of the line-of-sight veloci- 
ties. In Figures [3] and 2] we measure the moments of the 
line-of-sight velocities on inclined galaxies with an incli- 
nation 5 of i = 70° and then correct for the inclination 
by diving the velocities by sin(z). Maximal rotation ve- 
locities 7j rot are taken to be half the difference between 
the maximum and the minimum rotation velocity along 

5 IPizagno et~aLl pOOTl) measure the SDSS i-band Tully-Fisher 
relation, see Figure 15 1 on a sample of disk galaxies with a median 
inclination of i = 72° . 
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the slit. Central velocity dispersions er ccn are measured 
within a 1 kpc aperture. The measured value is increased 
by a factor 1.1 to account for the radial decline of the ve - 
locity dispersion (|Carollo et all 119931: iTreu et al.|[200lD . 
The quantity v/a refers to the ratio of maximum rota- 
tion velocity and central velocity dispersion. 

Cold gas is defined as the gas with a temperature of 
less then 3.2 x 10 4 K, unless noted otherwise. The neu- 
tral hydrogen mass is derived from the hydrogen mass 
fraction that is calculated within GASOLINE during the 
simulation. Our simulation does not include an explicit 
treatment of molecular hydrogen. Hence, in the analysis 
we do not distinguish between atomic and molecular hy- 
drogen components. Colors (all in the AB-system) are 
derived fr om stellar synthesis models , typically from the 
models of iBruzual fc Charlotl (|2003l ). but we ha v e also 
tested the single stellar populations of iMarastonl (2005) 
finding that the latter reddens the B — I color at z = 0.1 
of galaxies with B — I< 1.1 by up to 0.1 mag. Hence un- 
certainties in the stellar population models are unlikely 
to have strong effects on our results, but can nonetheless 
change colors and magnitudes by a few percent. 

We reduce the effect of star formation in the cen- 
tral region on the colors by fixing the mass-to-light ra- 
tio from within one softening length to its value at one 
softening length. In order to estimate the effect of a 
very young (< 100 Myr) stellar component, we have 
te sted a simplified versio n of the dust-obscuration model 
of ICharlot fe Falll (|2000t ). namely, reducing the flux of 
the very yo u ng st ellar particles to zero as done in, e.g., 
iNaab et all (|2007|) . but not applying any further dust 
corrections. We find that galaxies with B— I colors below 
~ 1 are typically reddened by 0.05-0.1 mag and redder 
galaxies are affected on a less than 0.05 mag level. Given 
that the corrections are overall rath er small we will re- 
fer to the single stellar populations of lBruzual fc Charlotl 
(2003) without a flux reduction of very young stars. 

The smoothed particle hydrodynamics technique al- 
lows us to follow the temperature history of individ- 
ual gas parcels. We track the temperature history of 
each gas particle over all snapshots in order to assess 
whether accreted gas has been cooled from a hot halo 
or whether it stayed cold since the beginning of the sim- 
ulation. We determine for each gas particle its maxi- 
mum temperature over its history, but ignore the times 
when a gas particle is less than 30 kpc away from the 
satellite center in order to avoid that its m aximum tem- 
peratu re is set by supernova feedback, see iBrooks et al.1 
(2009). We identify gas particles with a maximum tem- 
perature above (belo w) 2.5 x 10^ K as contr ibutors to 
hot (cold) accretion (jKeres et al.ll2005l [20091. Despite 
its simplicity this temperature-based split into hot/cold 
accretion gives simi lar results compare d to more sophis- 
ticated treatments ([Brooks et al.ll2009h . Since each star 
particle is spawned from a unique gas particle we can 
also use the maximum temperature of the spawning gas 
particle to asses whether a star particle has been formed 
from gas that was accreted hot or cold. 

To obtain the net, instantaneous accretion rates of gas 
and stars onto a satellite galaxy we first determine its 
center-of-mass velocity from all particles within a sphere 
of 20 kpc. Radial velocities v r of gas and star particles 
are then computed in the rest frame of the satellite. We 
then estimate the net mass flux of, e.g., gas with a max- 




FlG. 1. — Simulated group at z = 0.1. A massive central galaxy 
is surrounded by a large number of massive satellite galaxies. The 
image is 1.2 Mpc X 1.2 Mpc across and shows the projected mass 
densities of dark matter (blue; from 1.45 to 1450 Mq pc -2 ), neutral 
hydrogen (green) and stellar matter (yellow; both from 0.36 to 365 
Mq pc — 2 ). The white circle indicates the virial radius R200 = 396 
kpc of the group. The properties of the satellites with a stellar 
mass > 10 10 Mq (numbered) are shown in Figure [2] The galaxies 
are either entering the group for the first time (e.g., galaxy 1), or 
orbit within the virial radius of the group, or already crossed the 
group and are now at an apocenter (e.g., galaxy 4). 

imum temperature below 2.5 x 10 5 K, by summing over 
all particles of the relevant type in a spherical shell at 
r — 35 kpc distance from the center 

dM/dt ss Airr 2 p V" rrnvl / m l — 4irr 2 V" m,«[/V. 

i i i 

Here, p and V are the mean density and volume of the 
spherical shell. Its thickness is 10 kpc to ensure that 
it contains at least several hundred particles. For the 
accretion of bound material we only include the particles 
identified as bound to the satellite by AHF, otherwise we 
use all the particles in the shell. 

We quantify the mass loss of cold gas due to consump- 
tion and stripping by tracing forward in time the cold 
gas particles that are contained within 20 kpc at a snap- 
shot ~ 100-200 Myr before infall time. We consider that 
an initially cold gas particle may be kept as gas particle 
within a 20 kpc distance ( "kept" ) , that it may be kept af- 
ter being transformed into a star particle ( "consumed" ) , 
or that it leaves the 20 kpc radius ("stripped"). A small, 
but non-negligible fraction of the "kept" gas is heated, 
e.g., by supernova feedback. However, to first order the 
heating and cooling of gas should balance and hence the 
heated gas mass will be replaced by a similar gas mass of 
cooled gas. Hence, the "kept" gas mass should roughly 
correspond to the kept cold gas mass. The stated mass 
fractions and errors refer to the mean and the standard 
deviation of the various mass fractions for all z = 0.1 
gas-poor satellites. 

We measure single component Sersic indices n ser sic 
on face-on and edge-on mock images with GALFIT 
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Fig. 2.— Massive (M > 1O 1O M ) gal axies orbiting within the group at z = 0.1. Each composite image shows an cdgc-on view, i.e. 
a view along the intermediate axis of the reduced moment-of-inertia tensor of the stellar component within 10 kpc, of a satellite galaxy. 
The galaxies are arranged with respect to their color, Sersic index, rotational support, star formation rate, and gas fraction. Fundamental 
galaxy properties are shown at the top right of each panel (v/rr, B — I color, Sersic index, stellar mass in Mq, SFR in Mq yr — x , HI 
fraction). Galaxies 1, 2 and 6 are classified as normal disk galaxies, 3 — 5 as gas-poor disks and 7 — 13 as elliptical galaxies. A representative 
of each class is marked with a colored star. (Red star) This galaxy has a red color, hosts virtually no neutral gas, is almost spherical 
(sphericity of the stellar component c/a = 0.77), has a steep surface mass profile, and its kinematics is dominated by velocity dispersion. 
(Yellow star) The galaxy is a red, gas-poor, rotating disk galaxy (c/a = 0.38) with an /-band exponential scale length aj = 1.4 kpc. (Blue 
star) This galaxy hosts ~ 10 9 Mq of HI, has a bluer B — I color, and harbors a rotating stellar disk (c/a = 0.36) with aj = 2.3 kpc. Th e 
RGB color channels of the images correspond to the surface brightness in the rest- frame Bessel B, R and / filter bands ( Besselli 1 19901 ) ■ 
respectively, and range from 16 mag arcsec - 2 to 24.5 mag arcsec~ 2 . Green contour lines indicate the column densities of neutral hydrogen 
corresponding to 0.1, 1, 10, and 100 Mq pc -2 , while thin (thick) red contours correspond to column densities of 1 (10) Mq pc~ 2 of bound 
hot (T > 2.5 X 10 5 K) gas. The asymmetric gas contours in the two top left panels reveal the action of ram-pressure as the satellites move 
through the intragroup medium. The scale of all images is indicated in the bottom right panel. The green circle at the bottom-left corner 
of panel 13 has a radius of two gravitational softening lengths and indicates the resolution limit of the simulation. 



(jPeng et all |2002). In particular, we compute noise 
("sigma") maps based on the particle-per-pixel number 
in mock images of the surface mass density assuming a 
Poisson statistics. We exclude the central three soften- 
ing lengths from the fit and use the average value of the 
face-on and edge-on n scrs j c value. This approach sub- 
stant ially reduces the s purious effects of "overcooling" 
(e.g.. iMaver et al.l l2008) in the morphological classifica- 
tion, as it puts weight on the mass profile outside the 
very central regions. Our classification criterion, albeit 
coarse, allows us nonetheless to distinguish between the 
disk and elliptical satellites which are discussed in the 
main text. 

4. RESULTS 

4.1. A Hubble Sequence in a Galaxy Group 

The simulation was evolved down to redshift z = 
0.1. At this epoch, 13 well-resolved galaxies, with stel- 



lar masses > 10 10 M Q , orbit in, or near, the galaxy 
group, see Figure [TJ The stellar mass budget of the 
group is dominated by an early-type galaxy at its cen- 
ter, whose properties are discussed in detail elsewhere 
(|Feldmann et al.ll2010D . The mass range of these satellite 
galaxies straddles across the transition mass of ~ 3 x 10 10 
Mq below (above) which disk (early-type) galaxies are 
more numerous than galaxies with early-type (disk) mor- 
phology. Indeed, Figure [5] shows that, at z = 0.1, 
the simulated satellite galaxies span a broad range in 
morphological properties, from very flat disks to al- 
most spherical spheroids, as well as in rotational sup- 
port Kot/CTcen ~ 0.4 - 1.3), color (B - I = 0.6 - 1.3 
mag), neutral hydrogen fraction (/ gas ~ 0% — 5%), and 
SFRs. Properties range from passively evolving galaxies 
to moderately star-forming systems, with SFRs that are 
typical of similar mass z ~ disks, i.e., ~ 1 M Q yr _1 . 
Bulge to disk decompositions remain challenging even 
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2 2.5 3 3.5 4 4.5 1.7 1 1 1 1 1 1 1 

1.16 log l0 ( a [km/s])- 1.21 log l0 (R eff [kpc] ) + 0.55 " 18 " 19 " 20 M " 21 " 22 " 23 

Fig. 3. — Comparison between simulated and observed galaxy scaling relations. Simulated values are inferred from the z = 0.1 snapshot; 
blue (empty) and red (filled) symbols correspo nd to satellites tha t are classified as disks and ellipticals, respectively. Left: projection of 
the fundamental plane of early- type galaxies (Graves ct al. 2009) • The gray stripe indicates the location of real galaxies and red dots 
are our simulated objects with an early-type morphology by redshift z = 0.1. Right: the Tully-Fisher (TF) relation. Star symbols show 
maximum stellar linc-of-sight rotation velocities of the disk satellites; triangles show the maximum line-of-sight rotation veloc ities of neutral 
hydro gen for the two disk satellites with / gas > 3%. The solid line shows the i-band TF relation from the SDSS survey llPizagno et al.l 
2007) and the shaded area indicates its scatter. HI rotation velocities and i-band luminosities of the satellites arc consistent with the TF 
relation. We verified (not shown) that this result also holds for dust-corrected S-band luminosities (Vcrhciicn & Sancisi 2001; Vcrheiicn 
r2^1;[Beirfc de J ong 2001). Error bars in the top left corners of each figure include statistical errors for individual satellites (including 
Poisson noise in the galaxy representations, binning errors and fitting errors), and an estimate of the systematic errors in converting 
simulated data into observables. These systematic errors include contributions from the fixed aperture measurements used to measure 
the half-light radii, from estimating the luminosities from stellar synthesis models, from neglecting dust effects, as well as from differences 
between kinematic data extracted from moments and from Gauss-Hermitian fits. 



at the sub-kpc resolution of our simulation. However, 
based on standard morphological classification criteria 
that use global galactic properties such as the ratio of 
stream-to-dispersion velocities, and the index of single- 
component Sersic fits to the surface brightness profiles, 
we identify six group members as rotationally supported, 
disk galaxies (v rot /a cen > 1, n Scrsic < 2.5). These 
galaxies show a range of colors and SFRs that range 
from those of normal star-forming spirals to those of 
passive spirals/SO galaxies. Specifically, three of the 
disk galaxies have red colors (B — I ~ 1.2 — 1.3), con- 
tain no significant amount of neutral hydrogen (< 1%), 
and have stopped forming stars (top right in Figure [2|). 
The three other disk galaxies, instead, have bluer colors 
(B — I~ 0.6 — 0.9), significant neutral hydrogen fractions 
(> 3%), and SFRs of ~ 1 Mq yr" 1 (top left in Figure 
[3]). The remaining seven satellites have properties that 
are typical of early-type (hereafter "elliptical") galaxies, 
i.e., n Scrsic > 2.5 and v rot /(T ccn < 1. 

4.2. Scaling Relations of the Simulated Satellite 
Galaxies 

Real galaxy populations are observed to obey funda- 
mental s caling laws, such as the fundamental p lane (FP) 
relation (jDressler et al.lll987HBender et alll992D a nd the 
Tully-Fisher (TF) relation (|Tullv fc Fisher! Il977t ). The 
FP for elliptical galaxies shows that most such galaxies 
populate a thin two-dimensional surface in the three- 
dimensional space of characteristic stellar velocity a, 
half-light radius r c g , and surface brightness within the 



half-light radius, E ff. The FP reflects the virial relation 
between masses, sizes, and velocities; the observed tilt 
relative to the virial plane is interpreted as a dependence 
of the mass-to-light ratio with galaxy mass. Similarly, 
the TF relation connects the luminosity of disk galax- 
ies and their rotation velocities. Observations indicate 
L oc V^ t with a = 3 — 4 depending on wavelength and 
how dust extinction is modeled. 

We use these scaling laws to assess how similar, in their 
global properties, our simulated objects are to real galax- 
ies. For this purpose, we extract from the last simulation 
snapshot (z = 0.1) the follo wing properti es for the satel- 
lites in the Bessel B band (Besscll 1990) and the Sloa n 
Digital Sky Survey (SDSS) i band (|Fukugita et al.l[l996T) : 

(1) circular effective radii r e g, that include half the total 
flux within projected 20 kpc. The flux from within one 
baryonic softening length is derived from the included 
mass and the mass-to-light ratio at one softening length. 

(2) The surface brightness which is proportional to the 
flux divided by r^ ff . (3) The flux- weighted stellar line-of 
sight velocity dispersion and rotation velocity. (4) The 
line-of-sight rotation velocity of neutral hydrogen for the 
two disk satellites with a sufficient cold gas content at 
z = 0.1, measured as for the stellar component. 

In the left panel of Figure [3J we analyze the position of 
the simulated galaxies within the FP. Four of the seven 
early-type galaxies in our simulated sample have FP 
properties of real galaxies of similar type. The remain- 
ing three are slightly more compact. Numerical effects 
or missing physical ingredients are likely to play a role in 



G 



Feldmann, Carollo & Mayer 



200 



150 



100 



50- 



-50 



-100 



-150 



- 2 °-°15 



-10 





d [ kpc ] 



10 



15 



Fig. 4. — Rotation velocities and velocity dispersions of the 
two z = 0.1 disk satellites with extended gas disks as function 
of radius. The kinematic properties of galaxy 1(2), see Figure [2] 
are shown in blue (red). Filled (empty) star symbols correspond 
to stellar rotation velocities (stellar velocity dispersions) along the 
slit. Dark (light) shaded areas show HI rotation velocities (HI 
velocity dispersions). The width of the shaded regions indicates 
the errors in the HI kinematic quantities. Errors are computed for 
each radial bin separately via bootstrapping. 
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Fig. 5. — Kinematic properties of massive z = satellites and 
their progenitors. Shown is the ratio of maximum stellar line-of- 
sight rotation velocity and central velocity dispersion vs. stellar 
mass M s tar(^). Plotted are the 13 satellites with mass in excess of 
10 10 Af Q at z = 0.1. The redshift z is indicated by color: z = 0.1 
(black), z = 1 (red), and z = 2 (green); Circles (stars) indicate 
galaxies with a sphericity larger (smaller) than 0.45; filled (empty) 
symbols denote galaxies that are classified at z = 0.1 as ellipti- 
cal (disk) galaxies. At z ~ 1 — 2 most progenitors are rotation- 
supported galaxies. 



generating this effect ()Meza et al.l 120 03; Feld mann et al.l 
'20 10;. For instance, spurious angular momentum loss 
due to limited resolution may increase th e amount of 
gas t hat overcools in the central region ([Maver et al.l 
2008). The la ck of feedback from an active galactic nu- 
cleus (AGN) (iTeyssier et al.l 1201 If) and of galactic out- 
flows (|Davel2009t lOppenheimer et al.ll2010D may increase 
the fraction of gas that cools and is converted into stars 
in the central region of the galaxy. Finally, the low star 
formation threshol d may reduce the effic iency of the su- 
pernova feedback (|Governato et al.ll2010D . 

In the right panel of Figure [3] we show how the ro- 
tation velocities of the simulat ed disk galaxies com pare 
with the observed TF relation (jPizagno et al.ll2007t ). At 
z = 0.1 only two of our six disk galaxies have extended 
gas disks and thus allow us to measure rotation veloc- 
ities based on HI. The luminosities and rotation veloc- 
ities of these two galaxies are in good agreement with 
the observed TF relation. The figure also demonstrates 
that rotation velocities derived from the stellar compo- 
nent are not a good proxy for the rotation velocity of 
the gaseous component, because a significant fraction of 
the kinetic energy of stars is in unordered motions. This 
latter point can also be seen in Figure 01 where we show 
the radial dependence of the HI and stellar kinematics, 
respectively, of the two disk galaxies with extended gas 
disks. Clearly, neutral hydrogen has a larger rotation 
velocity and a smaller velocity dispersion than the stel- 
lar component. The figure further shows that one of the 
galaxies has a "flat" rotation curve both in HI and in the 
stellar component. The gas disk of the other galaxy is 
highly asymmetric due to ram-pressure stripping and we 
can measure its HI rotation velocity out to about ~ 5 
kpc. At larger radii we observe a declining rotation ve- 



TABLE 1 

The Classification of z = 0.1 Group Satellites with 
Stellar Masses above 10 10 Mq and Their Main Progenitors 



Satellite Classification 


z = 0.1 


z = 0.5 


2 


= 1 


2 = 2 


# Red 


9 


4 










# Blue 


4 


9 




13 


13 


# Gas-poor 


10 (9) 


5 (2) 


1 


(1) 


(0) 


# Intermediate-gas 


3 (3) 


8 (6) 


9 


(3) 


2 (1) 


# Gas-rich 


0(1) 


(5) 


4 


(9) 


11 (12) 


# Quiescent 


9 


4 




1 





# Star-forming 


4 


7 




6 


1 


# Strongly SF 





2 




6 


12 


# Quiescent (total) 


8 


1 










# Star-forming (total) 


4 


8 




4 





# Strongly SF (total) 


1 


4 




9 


13 



Note. — The rest frame B — I color is used to define red 
[B — I > 1) and blue (B — I < 1) galaxies. Galaxies are defined 
as gas-rich if they have a Hl-to-stellar mass ratio or (for the val- 
ues given in brackets) a cold gas-to-stellar mass ratio in excess of 
10%. Gas-poor galaxies have / gas <1%. Galaxies with / gas in the 
range 1%-10% are referred to as intermediate-gas systems. Galax- 
ies with specific SFRs (SSFRs) < 0.02 Gyr — 1 are defined as quies- 
cent galaxies, and otherwise as star forming galaxies. Galaxies are 
furthermore considered to be "strongly star forming" if their SSFR 
is larger than the inverse of the age of the universe, i.e., SSFR > 
0.08 Gyr -1 at z = 0.1. The classification in the last three rows 
differs from the one in the previous three rows due to the inclu- 
sion of the star formation occurring within the unresolved central 
softening length ei, ar . Note that at z ~ 1 — 2 most progenitors are 
blue, gas rich, star-forming galaxies. 
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Fig. 6. — Assembly histories since z = 2 of all z = 0.1 group 
members with a stellar mass above 10 10 Mq. Red solid (blue 
dashed) lines show the stellar mass of 2 = 0.1 elliptical (disk) 
satellites as function of time (bottom axis) and redshift (top axis). 
Filled (empty) symbols indicate the time when the progenitor of 
a 2 = 0.1 gas-poor (intermediate-gas) satellite enters the virial 
radius of the group. Sudden increases in stellar mass are caused 
by merger events, often before infall into the group, while close 
pericentric passages that lead to tidal stripping reveal themselves 
through a step- wise stellar mass loss occurring after their entrance 
in the group. The gas-poor disk satellites do not grow significant 
stellar mass after they enter the group, due to quenching of their 
star formation. One of the elliptical satellites does not show signs 
of significant merging since z ~ 2.2 (red triangle). This galaxy 
assembles most of its stellar mass early on, before such epoch, i.e., 
much earlier than most disk satellites, and, not surprisingly, this 
results in a massive, very compact spheroidal. 



locity in the stellar component. 

We conclude that, overall, the simulated galaxies 
match the fundamental scaling laws obeyed by the < 
L* z — galaxy population reasonably well. Of 
course, the match between real and simulated galax- 
ies is neither expected nor observed to be yet per- 
fect. Yet, we have come a long way in reproduc- 
ing the rich diversity of the z — galaxy population, 
from disks to ellipticals. This is achieved with nu- 
merical recipes for star formation and stellar feedback 
which were originally tuned to produce a Milky- Way - 
type galaxy (jStinson et al.ll2006t iGovernato et al(l2007l ). 
Despite this success in reproducing galaxies with ap- 
proximatively realistic properties, we should point out, 
however, that the stellar mass to halo mass 6 ratios of 
the group satellites are similar to the ones obtained in 
previous cosmolog ical simulations of g alactic halos (~ 
20% - 80%): e.g.. lAbadi et all (120031): IGovernato et all 
(|2007h : IScannapieco et all (|2009f l: lAgertz et all (|2011f )7 
Halo abundance matching approaches based on the ob- 
served stellar mass function typically predict stellar mass 
to halo mass ratios that are lower by a factor of a 
few, e.g., iGuo et all (|2010f ). This mismatch may have 
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6 The latter is defined as the mass of the halo when the satellite 
first enters the virial radius of the group. 



7. — Events in the life of the group satellites. The x-axis 
indicates the time of the first pericentric passage i pcr i. The y- 
axis indicates the following events: the time when the satellite first 
crosses the virial radius of the group (shown using blue squares), 
the time when the Hl-to-stellar mass ratio drops below 1% (green 
triangles) , and the time when the B — I color increases above 1 (red 
circles). Note that some satellites have not turned into gas-poor 
or red galaxies by 2 = 0.1 and, therefore, for these galaxies green 
triangles or red circles are missing. Galaxies that are identified as 
disks (ellipticals) at 2 = 0.1 are shown with empty (filled) sym- 
bols. Linear fits with unity slope are shown as colored solid lines, 
and the fit errors as dotted lines. The satellites approach their 
first pericentric passage typically 0.4 ±0.1 Gyr after infall into the 
group and, by then, turn into gas-poor (/gas < 1%) galaxies. On 
average 0.6 ± 0.5 Gyr later the galaxies become red. The error of 
a typical data point that arises from the limited time resolution of 
the simulation is indicated at the top left. The arrow at the top 
right highlights a gas-rich disk that has not yet completed its first 
pericentric passage. A Spearman rank test shows a correlation, at 
the level of 2.8ct and 2.2<r, respectively, between pericentric time 
and both the time when / gas drops below 1% and the time when 
the galaxy turns red. 

its origin in an over-efficient conversion from gas to 
stars in simulations indicating that the feedback pre- 
scription is incomplete and that additional mechanisms 
that remove gas from the galaxies, e.g., radiation pres- 
sure (Hopkins et al. 2011), or supernova driven outflows 
(jGuedes et alJl20ill ). are required. 

4.3. The Drivers of the Transformation Processes 

At z > 1, the most massive progenitors of the z = 
group satellite galaxies are blue (B — I < 1), strongly 
star- forming (specific SFR > tjj 1 ), gas-rich (/ gas ~ 20%), 
and rotationally-supported (v vo t/a CC n > 1), see Table Q] 
and Figure [SJ In the simulation we can trace when, 
and how, the progenitor disks of the passively evolving 
spheroids and "gas-starved" disks, exhausted and/or ex- 
pelled their fuel for star formation, and when and how 
they transformed their morphologies (in those cases they 
did). 

Major galaxy mergers are in most cases responsible 
for the morphological transformation of disk galaxies 
into systems with an elliptical morphology. Remark- 
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Fig. 8. — Gas accretion history of (left) a z = 0.1 gas-poor disk satellite (3 in Figure 2), (middle) a z = 0.1 intermediate-gas disk satellite 
(2 in Figure 2), and (right) a z = 0.1 elliptical satellite (13 in Figure 2). Shown are the instantaneous accretion rates of gas that was 
never heated above T t i lrcs = 2.5 X 10 K (cold accretion; blue line), and of gas heated above that threshold (hot accretion; red line). The 
dotted vertical line denotes the infall time into the group, and arrows indicate subsequent pericentric passages. The solid line includes only 
particles bound to the respective satellite, while the dashed line includes all particles (only shown before infall). The right panel also shows 
the mass accretion rates of stellar particles that are formed by gas particles once heated (cyan) or not heated (magenta) above Tthres- 
The gas accretion of disk satellite progenitors (left and middle panel) proceeds in a smooth fashion, but is very irregular for the 2 = 
elliptical satellite. The gas accretion of the gas-poor disk is already quenched at t ~ 7 Gyr, when it resides in a small neighboring group. 
The gas accretion of the intermediate-gas disk is terminated after infall to the main group. Just before the infall the density of the hot 
gas increases as the satellite enters the group halo, which is filled with ho t, shock- heated intr agroup gas. Shortly afterward, ram pressure 
stripping removes a large fraction of the satellite's own hot halo gas (McCarthy ct al. 2008). In the case of the elliptical satellite cold 
accretion dominates at early times (t < 4 Gyr). At later times, but before the fall into the group, the gas accretion history is dominated 
by mergers. At the latest epochs the group environment quenches any further gas accretion onto the satellite. 
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Fig. 9. — Fate of the cold gas in (left) a z = 0.1 gas-poor disk satellite (3 in Figure 2), (middle) a z = 0.1 intermediate-gas disk satellite 
(2 in Figure 2), and (right) a z = 0.1 elliptical satellite (13 in Figure 2). The cold gas particles contained within a sphere of 20 kpc radius 
around the satellite are identified in a snapshot before infall time (vertical dotted line) and traced forward in time. The solid black line 
shows the total mass of all successor particles of the selected cold gas particles. This mass is not conserved but slightly increases due to 
stellar winds from the satellite's own stellar population. The mass that remains within (< 20 kpc) the satellite is shown as dashed black 
line, while the mass of stripped (> 20 kpc) particles is shown with a dot-dashed black line. The other lines correspond to the amount 
of gas that remains and stays cold (dashed green line), that is transformed into stars that remain within the satellite (dashed cyan line), 
that is stripped (dot-dashed black line), that is stripped and heated (dot-dashed red line), that is stripped and stays cold (dot-dashed blue 
line), or that is transformed into stripped stars (dot-dashed cyan line). Most of the cold gas reservoir present before infall into the group 
is converted into stars within 1 Gyr or less, while a small fraction is stripped between infall time and the first pericentric passage (vertical 
arrow) . 



ably, and rather unexpectedly, however, these transfor- 
mations happen before galaxies enter the virial radius 
of the group. They produce six out of the seven ellipti- 
cals that are present in the group by the present epoch. 
Therefore, while groups have always been considered as 
the natural site for major galaxy mergers given their 
low velocity dispersion, we find that elliptical galaxies in 
groups are the result of pre-processing at z > 1, caused 
by interactions occurring when galaxies are still distinct 
objects embedded in their individual halos. The galaxies 



are also still gas- rich when they undergo such mergers, 
and they will lose their gas only later when they will fall 
into the virialized region of the group. Figure [6] shows 
that the progenitors of z — 0.1 elliptical satellites all suf- 
fered from (moderately dissipative, / gas ~ 10% — 20%) 
major (> 1 : 6) merger events over the z ~ 1 — 2 pe- 
riod. In contrast, the progenitors of z — disk satellites 
have a rather quiet and uneventful history: none of them 
experiences a significant (> 1 : 10) stellar merger after 
rcdshift z ~ 2. These galaxies grow their stellar mass by 
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Fig. 10. — Effects of pericenter passages in the evolution of the 
satellites. The figure shows the average cold gas transformation 
rates obtained after stacking all satellites together on normalized 
axes. The time axis is normalized in units of t pcr i — ii n f a ii> while 
the rates of each individual satellite are normalized to its maxi- 
mum value. Thick (thin) lines denote the median (lower/upper 
quartile) of the satellite population. Most of the cold gas that 
resides in satellites just before they infall into the group (green 
dashed lines; shown with the negative of the mass loss rate for 
better visualization) is either continuously transformed into stars 
(blue dot-dashed lines) or stripped (black solid lines) in pericen- 
tric passages. Some satellites show enhanced star formation rates 
close to pericenter time. However, overall, the conversion of cold 
gas into stars is strongly reduced after the first pericentric passage. 
Within a timescale of about i pcr i — Onfall after pericenter passage, 
the cold gas reservoir is exhausted and the mass loss rates become 
negligible. 

in situ star formation, which typically proceeds at a rate 
~ 1 Mq yr _1 at z < 1. The mergers that produce the 
elliptical morphologies naturally lead to a faster growth 
of stellar mass relative to in situ star formation. 

The group potential plays on the other hand a cru- 
cial role in removing the gas reservoir of all galaxies and 
truncating their star formation (Figure [7]). As a result, 
the photometric evolution of ellipticals into the passively 
evolving early-type galaxy population observed at the 
present epoch lags behind by several Gyr the time at 
which their spheroidal morphology is established in the 
mergers at earlier epochs. Likewise, the transformation 
of normal star-forming disks into gas-poor, "starved" 
red disks appears to be mostly a consequence of their 
infall into group environment. Indeed, the three disk 
galaxies that are gas-poor and quiescent by z — en- 
ter the group at earlier epochs than the three z = 
star-forming disks (respectively, z — 0.35, 0.49, 0.80 and 
z = 0.13,0.19,0.26), and are exposed for longer periods 
to the physical processes that are most effective in the 
group environment. Such processes are, precisely, sup- 
pression of gas accretion as well as ram pressure stripping 
of both the hot and cold gas reservoir; in situ star forma- 
tion is responsible for consuming a substantial fraction 
of the residual gas. In detail, a galaxy that plunges into 
the depth of the group potential is first observed to stop 




Fig. 11. — Formation of a z = 0.1 elliptical galaxy that lives 
today in the relatively dense potential of a ~ 10 13 Mq galaxy group. 
(Top left) Edge-on view of the z = 1.1 main progenitor of the 
elliptical galaxy 13 in Figure [2] the progenitor has both a stellar 
and a gas disk. (Top right) A gas-rich merger occurs ~ 400 Myr 
later, destroying the stellar disk. (Bottom left) Another 300 Myr 
later the stellar component of the galaxy has relaxed to an elliptical 
morphology while still showing blue colors. Also, a small gas disk 
begins to grow again, out of the remaining cold gas. (Bottom right) 
At z = 0.38 the galaxy completes its first pericentric passage. By 
this time the galaxy has lost most of its cold gas, its star formation 
is quenched, and the colors of the galaxy are now typical of "red 
sequence" , passively evolving galaxies. Column densities of neutral 
hydrogen are indicated by green contour lines which correspond to 
0.1, 1, 10, and 100 M pc -2 , respectively. 



accreting gas after the infall. Subsequently, a substantial 
amount of its hot halo gas is removed by ram-pressure 
stripping on its first orbit within the group, namely on a 
timescal e of a Gyr or so, consistent with a "starvation " 
picture (jLarson et all 119801: iKawata fc Mulchaevl I2008D . 
see Figure H Most of the cold gas (~ 90 ± 10%) is then 
lost within the following ~ 1 Gyr, through star formation 
(~ 70%) as well as nearly instantaneous gas stripping at 
pericenter passages (Figures GO and [TU]). For galaxies 
falling into the final virialized group halo, the first peri- 
center passage is reached 0.4±0.1 Gyr after infall into the 
group, i.e., the characteristic crossing time of the group 
(Figure El). 

In Figure [TT] we illustrate, as an example, how the 
transformation into a gas-poor, quiescent, red, velocity- 
dispersion supported elliptical galaxy takes place as a 
two-step process, starting from gas-rich, star-forming, 
blue, rotating galaxies. First, at z ~ 1.1 and at ~ 1 
Mpc distance from the group center, a gas-rich merger 
with a stellar mass ratio ~ 1 : 2 destroys the stellar disk 
of the more massive progenitor. The remnant relaxes in a 
few hundred Myr to an elliptical morphology, but its col- 
ors remain blue and it still hosts a substantial amount 
(~ 10 9 Mq) of HI within 20 kpc. During the merger 
the S FR is substantially enhanced (Mih os fc Hernquistl 
1996); at its peak, it is about a factor five higher than the 
average rate that the galaxy experiences throughout its 
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Fig. 12. — Evolution of morphological, kinematic, and photometric properties of (left) a z = 0.1 gas-poor disk satellite (3 in Figure 2), 
(middle) a z = 0.1 intermediate-gas disk satellite (2 in Figure 2), and (right) a z = 0.1 elliptical satellite (13 in Figure 2). Shown are 
the line-of-sight velocity dispersion (magenta line with error bars), the maximum of the line-of-sight rotation velocity (red line with error 
bars), v/a (blue line with error bars), the non-sphericity 1- c/a (black solid line), the amount of cold gas (< 3.2 X 10 4 K) within 20 kpc 
(green dot-dashed line), the star formation rate within 10 kpc (upper red dotted line), within 10 kpc excluding one softening length (red 
dot-dashed line), within 10 kpc excluding two softening lengths (lower red dotted line), and the B — I color (magenta shaded area; the 
upper and lower limits are measured within 10 kpc excluding and including the central two softening lengths, while the magenta line is 
measured by excluding the central softening length). The red arrow at the top left is the time when half of the stellar mass is formed. The 
vertical dotted line denotes the time of infall into the group and the black arrows on the bottom right denote pericentric passages. The 
vertical shaded areas in the right panel denote prominent mergers (red) or a fly-by (blue). 



history. Thereafter the HI reservoir and the SFR gradu- 
ally decrease. After the first pericenter, the SFR and the 
HI content become negligible and the galaxy turns into 
a red, passively evolving el liptical (~ L* in th e SDSS 
z = 0.1 redshifted r-filter; iBlanton et al.l 120031 ) with a 
final B - I = 1.34. 

We note that one of the seven z = elliptical satellites 
does not undergo any substantial merger since z ~ 2.2 
(red triangle in figure Figure EJ) . This galaxy assembles 
most of its stellar mass early on, i.e., much earlier than 
most disk satellites. Not surprisingly, this results in a 
massive and very compact spheroidal system. We specu- 
late that the very diverse stellar density and velocity dis- 
persion pro perties observed in the z ~ 2 passive galaxy 
population (jvan Dokkum et all [20091 could be, at least 
in part, explained by a spread in the assembly epochs of 
its constituent galaxies. 

4.4. Environmental Effects before Infall into the Group 

Although the ultimate quenching of galaxies is caused 
by stripping of the gas reservoir and the resulting sup- 
pression of star formation after galaxies enter the group 
potential, we note that their star formation histories start 
being affected by environmental effects also prior to their 
infall. For most galaxies, the SFR peaks around z ~ 2 
at a distance of ~ 1 Mpc (~ turn around radius) from 
the center of the assembling group. Afterward, a differ- 
ent behavior is seen in the progenitors of the present- 
day population of, respectively, disk and elliptical galax- 
ies. In the progenitors of disks, star formation declines 
smoothly from high to low redshift until truncation en- 
sues, reflecting the fact that these galaxies evolve in near 
isolation before joining the group. In contrast, the SFR 
of the progenitors of ellipticals, which are gas-rich, mas- 
sive disk-like objects, shows high variability, with bursts 
caused by galaxy interactions and mergers (Figure [T2|). 

As similar dichotomy occurs for gas accretion, which 
happens predominantly in the "cold mode" for all galax- 



ies prior to infall onto the group. In galaxies which re- 
main disks down to z = 0, the accretion rate decr eases 
gradually with time (jSommer-Larsen et al.l 120031 ). see 
Figure |U In the early- type galaxies, the rate is instead 
rather bursty and strongly modulated by the mergers 
these systems endure. In some galaxies, the decrease in 
the cold accretion rate is exacerbated by "truncation" of 
the streams due to the reduction of the tidal radius of 
the satelli tes (driven by the tidal effects of the group halo 
potential: lHahn et al.l 12009). This leads to quenching of 
gas accretion already at a few virial radii away from the 
group center. In at least one of the gas-poor disks (Fig- 
ures [H and IT!?]) , the supply of cold gas is already termi- 
nated when it resides in the vicinity of a ~ 5 x 10 11 M© 
halo that falls into the main group at a substantially later 
time. 

5. CONCLUSIONS 

Our work shows that the formation of < L* elliptical 
galaxies and "starved" red disk galaxies in high-density 
environments such as galaxy groups is the result of the 
combined effect of merging, gas removal and gas con- 
sumption processes as well as halted gas accretion natu- 
rally occurring as the potential of the group is assembled 
hierarchically in a ACDM universe. Our simulations do 
not include feedback from AGNs and thus show that it 
is not an essential ingredient in order to form elliptical 
galaxies with quenched star formation in z ~ galaxy 
groups. 

Specifically, we find that elliptical galaxies are formed 
by mergers at high redshift prior to the formation of the 
final group, which is consistent with the fact that ellipti- 
cal galaxies are found to be already in place at high red- 
shift e ven in the field (ICimatti et al.ll2006l ; IScarlata et al.l 
120071: iBolzonella et al.l 120101) . The lag of photomet- 
ric transformations relative to the morphological trans- 
formations seen in our simulation indicates that many 
of the "excess" blue L* galaxies observed at z ~ 1, 
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which are thought to become L* ellipticals, but which 
appear to have an irregular morphology at blue rest- 
frame wavelengths (Scarlata ct al. 20 07tlBolzonella et al.l 
2010), should show an elliptical morphology when ob- 
served at longer rest-frame wavelengths. This is a 
straight prediction of our simulations since at the time 
when their old stars acquire an elliptical morphology 
through mergers, the simulated galaxies still host a star- 
forming, gaseous component. This also indicates that 
the tight correlation of colors and Hubble types, as well 
as the standard classification based on the Hubble se- 
quence, does not hold for galaxies at redshifts higher 
than about 1.5 — 2, in good agreement with state-of-the- 
art studies of galaxy m orphologies at such high redshifts 
(|Cameron et al.ll2010[) . 

In addition, our model predicts that the star formation 
and gas accretion rates of progenitors of elliptical galax- 
ies in the considered mass range (~3x 10 10 M Q ) should 



be highly irregular due to gas-rich galaxy interactions 
and mergers at intermediate redshifts z > 1. Therefore, 
these galaxies should show a large scatter in the ages of 
their stellar populati ons and differ from s i ngle, high red- 
shift burst models (jGallazzi et al.l l2005t iThomas et al.l 
2005). Future multi- wavelength observations capable of 
characterizing with unprecedented detail both the SFR 
evolution and the evolution of the gas content of galaxies, 
atomic and molecular, as a function of redshift, such as 
the Atacama Large Millimeter Array, the James Webb 
Space Telescope, and the Square Kilometer Array, will 
be able to test this scenario. 
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